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The prospects for indirect detection of dark matter at tlie Galactic center with gamma-ray exper- 
iments like the space telescope GLAST, and Air Cherenkov Telescopes like HESS, CANGAROO, 
MAGIC and VERITAS, depend sensitively on the mass profile within the inner parsec. We calculate 
the distribution of dark matter on sub-parsec scales by integrating the time- dependent Fokker-Planck 
equation, including the effects of self-annihilations, scattering of dark matter particles by stars, and 
capture in the supermassive black hole. We consider a variety of initial dark matter distributions, 
including models with very high densities ("spikes") near the black hole, and models with "adia- 
batic compression" of the baryons. The annihilation signal after 10^" yr is found to be substantially 
reduced from its initial value, but in dark matter models with an initial spike, order-of-magnitude 
enhancements can persist compared with the rate in spike-free models. 



PACS numbers: 

There is compelling evidence that the matter density of 
the Universe is dominated by some sort of non-baryonic, 
"dark" , matter, the best candidates being weakly inter- 
acting massive particles 0, Q • Numerical iV-body sim- 
ulations suggest dark matter (DM) density profiles fol- 
lowing broken power laws, p oc r~'^ , with 7 « 3 in the 
outer parts of halos and 1 < 7 < 1.5 ("cusps") inside the 
Solar circle. Although these profiles reproduce with suf- 
ficently good accuracy the observed properties of galac- 
tic halos on large scales, as inferred by rotation curves, 
little is known about the DM distribution on smaller 
scales, where the gravitational potential is dominated by 
baryons. The situation at the Galactic center (GC) is fur- 
ther complicated by the presence of a supermassive black 
hole (SBH), with mass ~ 10^-^ M© 0|, whose sphere of 
gravitational influence extends out to ^ 1 pc. 

The prospects for indirect detection depend crucially 
on the distribution of DM within this small region. The 
flux of gamma-rays from the GC, from the annihilation 
of DM particles of mass m and annihilation cross section 
in the non-relativistic limit av, can be written: 

<,.m,E)^<^/-^(j^) fi^)'j.oA17 (1) 

dE \{av)thj \ m J 

where $0 — 5.6 x 10^^^cm~^s~^ and and {av)th ~ 
3 X 10~^^cm"^s~^ is the value of the thermally averaged 
cross section at decoupling that reproduces the observed 
cosmological abundance of Dark Matter (although in 
presence of resonance effects like co-annihilations, the 
correct relic abundance can be achieved with smaller 
cross sections). For more details and a review on DM 
candidates and detection see e.g. Refs.0, ^An is 
a factor containing all the information on the DM pro- 
file 0: 

Jao = K{An)-^ [ di, [ p^dl, (2) 
J An J4> 

where dl is the distance element along the line of sight 
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at angle ip with respect to the GC, Ai7 is the solid 
angle of the detector, and K is a. normalizing factor, 
K-^ = (8.5kpc)(q.3GeV/cm3)2. We denote by Jg and 
J3 the values of J when Afi = 10~^sr and 10~^sr re- 
spectively; the former is the approximate field of view of 
the detectors in GLAST Q and in atmospheric Cerenkov 
telescopes like VERITAS and HESS 0, while the 
larger angle corresponds approximately to EGRET |^. 
DM densities that rise more steeply than p oc r~'^/^ near 
the GC imply formally divergent values of J, hence the 
predicted fiux of annihilation products can depend sen- 
sitively on any physical processes that modify the DM 
density on subparsec scales. Although the analysis of DM 
indirect detection is usually performed under simplifying 
assumptions on the DM profile - extrapolating the re- 
sults of numerical simulations with power-laws down to 
subparsec scales - several dynamical processes may in- 
fluence the distribution of DM at the GC, including the 
gravitational force from the SBH jO| , adiabatic compres- 
sion of baryons [3| , and heating of the DM by stars |0 . 

Here, we focus on the evolution of the annihilation sig- 
nal due to two physical processes that are almost certain 
to strongly influence the form of the DM density profile 
near the GC: DM self-annihilations; and gravitational 
interactions between DM particles and stars. Both pro- 
cesses act on a similar time scale (^ 10^ yr) to modify 
p{r) on the sub-parsec scales that are most relevant to 
the indirect detection problem. While these two mecha- 
nisms both tend to lower the DM density, we find that 
interestingly high densities can persist over a particu- 
lar range of (m, av) values. The time-dependent profiles 
discussed here may also have important consequences for 
the prospects of observing an extra-galactic gamma-ray 
background. 

Let /(r,v,i) be the mass density of DM particles in 
phase space and p{r, t) their configuration-space density, 
with r the distance from the GC, i.e. the distance from 
the SBH. We assume an isotropic velocity distribution, 
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/(r, V, t) = f{E, t), where E = -v'^/2 + (f){r) is the bind- 
ing energy per unit mass and —(j){r) is the gravitational 
potential, which includes contributions from the stars in 
the Galactic bulge and from the SBH. We assume that 
(f) is fixed in time, i.e. that the mass of the SBH has not 
changed since the epoch of cusp formation, and that the 
stellar distribution has also not evolved. The first as- 
sumption is commonly made based on the observed, very 
early formation of massive black holes (e.g. The 
second assumption is motivated by the expectation that 
the stars should reach a coUisional steady state around 
the SBH, the so-called "Bahcall-Wolf" solution 13], in 
a time of ~ 10^ yr. The observed stellar distribution at 
the Galactic center, p^, ~ ^-i-* jg slightly shallower 
than the Bahcall-Wolf solution for a single population 
but is generally believed to be consistent with a coUi- 
sional steady state given uncertainties about the stellar 
mass spectrum We accordingly set cx r~^-^ 

and fix its normalization to match the observed stellar 
density at ~ 1 pc from the SBH The stellar phase- 
space mass density /* {E) is then uniquel y d etermined by 
pj,(r) and 0(r) via Eddington's formula jl6lll7l |. 

We describe the evolution of / via the orbit-av erag ed 
Fokker-Planck equation including loss terms (e.g. [18]): 
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~f{E,t)ME), 
Dee{E) = Ui^^G'^m^lnK x 

q{E) / dE'ME') + / dE'q{E')U{E') 



(3a) 
(3b) 



Here p{E) = 4\/2 JJ— drr'^^4>{r) - E is the phase 

space volume accessible per unit of energy, p{E) = 

—dq/dE, and In A is the Coulomb logarithm |l8l |. 
We have assumed that the spectrum of stellar masses 
n{fni,)dm-i, is independent of radius; then m^, — 

The first term on the right-hand side of Eq. 3a de- 
scribes the diffusion of DM particles in energy space due 
to heating via gravitational encounters with stars. Near 
the SBH, the characteristic heating time is nearly in- 
dependent of energy and radius and is given approxi- 
mately by Thcat = 1-25 X lO^yrAf.^^^^rl^/^Tn^"^ with 
M, = M,/(3 X IO^Mq); fh = r/,/(2pc), where rh is the 
"gravitational influence radius" of the SBH, defined as 
the radius of the sphere containing a mass in stars equal 
to twice M,; and rfii, = rni,/MQ In what follows 

we set M, = fh = rfii, — 1 and define r = t/Thcat; the 
age of the Galactic bulge, ~ 10 Gyr, then corresponds to 
T « 10. (The most recent estimates of M, are slightly 
highe r .IQJ ; we adopt M, = 1 for consistency with earlier 
work 

The collision term Vcoii has two potential contributors: 
self-annihilations, and interaction of DM particles with 



baryons. The self-annihilation term is given locally by 
ly = m^^ pav. The orbit-averaged rate i/coii that appears 
in Eq. 3a is 



i^coll(-B) = 



J h'r'^v{r, E)dr 
J r^v{r, E)dr 



(4) 



where v{r, E) — {(j>{r) — E) and the integrals are from 
to rinax(£^)- Expressing p in terms of /, we can write 
the orbit-averaged self-annihilation term as 



p{E)vco\\{E) = Z2'K{(jv)m-^ x 



nE nOO 

/ dE'f{E')C{E,E') / dE'f{E')C{E',E) 
Jo Je 

C{E,E'), 



(5a) 



drr^ y/(t){r) - Ey/(j){r) - E{5h) 



Self-annihilations limit the DM density roughly to p « 
m/{av)t Hill; for m = 50 GeV, av = lO'^^ cm^s^^ 
and t = 10 Gyr, p <2x 10^ Mqpc^^ , which would imply 
that self-annihilations are important at r < 10^'^ pc if 
p - and p(r/,) = IOOMq pc'^. 

The final loss term in Eq. 3 represents scattering of 
DM particles into the SBH This term, which we 

include, is important at radii r < r^- The loss rate 
varies only logarithmically with the SBH's capture ra- 
dius, which we set to 2GMt/c^. 

In what follows, we assume rhi, = IMq, consistent with 
our limited knowledge of the stellar mass spectrum near 
the Galactic center [ll| . We note that both the first and 
third terms on the right hand side of eq. 3a depend in 
the same way on m*; thus, varying m^, has the effect of 
changing the relative time scales for DM-star scattering 
and self-annihilations. Since changing av has the same 
effect, we do not vary in what follows. 

Eq. 3 was advanced in time via a backward differenti- 
ation scheme coupled with the method of lines to reduce 
the partial differential equation to a system of ODEs 
A variable time step was employed, such that the frac- 
tional change in / in one time step was less than 1% at 
every value of E. For very high values of av/m, the ini- 
tial / was truncated such that the annihilation time was 
never shorter than ~ 10^ yr. 

We adopted a wide range of initial conditions for the 
DM distribution (Table I). Baryon-free simulations of 
DM clustering suggest a power-law distribution in the 
inner parts of galaxies, p (x r~'^'=, a "cusp," with 7^ « 1 
(these models are labelled "N" in Table I). The most 
recent simulations (see also Refs. jHli^]) suggest 

a power-law index that varies slowly with radius, but the 
normalization and slope of these models at r « r/i are es- 
sentially identical to those of models with an unbroken, 
p cx power law inward of the Sun. We took Rq — 8.0 
kpc for the radius of the Solar circle . 

Since the total mass budget of the inner Galaxy is 
dominated by baryons, the DM distribution is likely to 
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TABLE I; Properties of the halo models. "N" and "A" stand 
for NFW and adiabatically contracted profiles, respectively. 
The subscripts "c", "sp" are for profiles with core and spike 
respectively. Core radius is in units of ru ■ Density at Rq is 
in units of GeV cm~'^. J3 and J5 are values of J averaged over 
windows of solid angle 10^^ sr and 10^^ sr respectively and 
normalized as described in the text. The final two columns 
give J in evolved models for av = Q (no annihilations), and for 
{av^m) = (3 X 10~^^cm^s~^, 50GeV) (maximal annihilation 
rate), respectively. 















logio (J5) 






7c 


-lap 


rc 


p[Rq) 


r = 


r = 10 


r = 10 


TV 


1.0 






0.3 


2.56 (3.51) 


2.56 (3.50) 


2.56 (3.50) 




1.0 




10 


0.3 


2.54 (3.33) 


2.54 (3.33) 


2.54 (3.33) 




1.0 


2.33 




0.3 


9.21 (11.2) 


3.86 (5.84) 


2.56 (3.52) 




1.0 


2.29 


10 


0.3 


6.98 (8.98) 


2.61 (3.88) 


2.54 (3.33) 


A 


1.5 






0.5 


5.80 (7.75) 


5.26 (7.03) 


5.23 (6.98) 


A, 


1.5 




10 


0.5 


4.96 (6.27) 


4.96 (6.27) 


4.96 (6.27) 




1.5 


2.40 




0.5 


14.8 (16.8) 


9.25 (11.3) 


5.25 (7.02) 




1.5 


2.29 


10 


0.5 


9.99 (12.0) 


5.21 (6.96) 


4.96 (6.27) 



have been influenced by the baryonic potential and its 
changes over time. The "adiabatic-growth" model ^3 
posits that the baryons contracted quasi-statically and 
symmetrically within the pre-existing DM halo, pulling 
in the DM and increasing its density. When applied to a 
DM halo with initial 7^ ~ 1, the result is a halo profile 
wit h « 1.5 inward of Rq and an increased density at 
Rq 123,113,1^. Adiabatically contracted halo models are 
labelled "A" in Table I. Alternatively, strong departures 
from spherical symmetry during galaxy formation might 
have resulted in a lower central DM density. For instance, 
the DM density following a merger is a weak power law, 
o_~ r-Tin^ ~ 0.5, inside a radius rc ~ 10 — lOOr^ 
[23, 113 • Models with the subscript "c" in Table I have 
p cx r~^/'^ inside a radius Tc = lOpc. 

We also considered modified versions of each of these 
DM profiles that included a density "spike" around the 
SBH; these models are denoted by the subscript "sp" in 
Table L The inner DM density in the spike models fol- 
lows the steeper power law that would result from grad- 
ual growth of the SBH to its current mass at a fixed 
location. We set p — p{rb){r/ri,)^'>'='p for r < ri, with 
jsp = 2-1-1/(4 — 7) ^-iid 7 the power-law index of the 
core or cusp, and r?, = 0.2r/j |9|,|3J|. It is unclear whether 
spikes can survive at the centers of all galactic halos, 
since dynamical effects such as off-center formation of 
the SBH and binary black hole mergers would tend to de- 
stroy high density regions [i^,!!^!- However, the Milky 
Way is unlikely to have experienced a "major merger" 
(a merger with another galaxy of comparable mass ) in 
the last 10 Gyr, and the existence of a stellar "cusp" [lj| 
further strengthens the case for a dark matter spike at 
the Galactic center 'X{ . 



In order to evaluate the influence of annihilations on 
the evolution of the DM profile, we first investigated two 
extreme cases in the framework of typical DM candidates 
like neutralinos or Kaluza-Klein particles 01 . In the first 
extreme case, in order to maximize the ratio av/m, we 
assumed a cross section av ~ {av)th = 3 x 10^^^ cm^ 
s~^ and a mass of 50 GeV. Higher values of the annihi- 
lation cross section, though possible, would imply a low 
relic density, making the candidate a subdominant com- 
ponent of the DM in the Universe, a case we are not 
interested in here. The lower limit on the mass strictly 
applies only to neutralinos in theories with gaueino and 
sfermion mass unification at the GUT scale |35j, while 
the limit on the mass of KK particles is higher. The 
second extreme case assumes no annihilations, as in the 
limit of very small cross sections, or very heavy particles. 
Table I gives values of J3 and J 5 at t = 10 for each 
of our DM models and for both of the extreme particle 
physics models. The J-values depend appreciably on the 
particle physics model only when the initial DM density 
has a "spike" around the SBH; in other cases the central 
density is too low for annihilations to affect J. Partic- 
ularly in the case of maximal crw, the final J values are 
found to be modest, logj^g '^3 ^ ^-^ and logj^p J5 < 7.0, 
compared with the much larger values at r = in the 
presence of spikes. 

We also carried out integrations for the set of bench- 
mark models derived in [36] in the framework of min- 
imal supergravity (mSUGRA). Although other scenar- 
ios (supersymmetric or not) predict different parameters 
for the DM candidate, the values of av/m are often ap- 
proximately the same. Light DM candidates [l^l, for 
example, have masses smaller than 20 MeV if they are to 
be responsible for the 511 keV emission from the Galac- 
tic bulge ^ , but they also typically have cross sections 
much smaller than the thermal cross section in the early 
universe, which implies that av/m falls again in the same 
range discussed above. Heavy candidates, like those pro- 
posed to explain the HESS data [13^3' have masses in 
the 10-20 TeV range, and thermal cross sections, so that 
they fall again in the same range of av/m. Fig. 1 shows 
the final DM density profile for each of the benchmark 
models, starting from DM models Nsp and Asp] the lat- 
ter model is the "adiabatically contracted" version of the 
former. While adopting the maximal annihilation rate 
effectively destroys the spike and produces J values as 
low as those of spike-free models, other benchmark mod- 
els with smaller av/m result in strong enhancements in 
J. 

Fig. 121 shows the evolution of the dark matter den- 
sity at a radius of lO^^r/i « 2 x lO^^pc, starting from 
a p ~ j,-2.33 gpjjjg (p ^ ^-1 cusp). Two values were 
taken for the initial density normalization at r = r/j, 
— (10, lOO)Af0pc~^. The self-annihilation term 
in Eq. (|Sa|l was computed assuming m — 200 GeV, 
av = lO^'^^cm'^s"^. The early evolution is dominated 
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FIG. 1: Evolved DM density profiles at r = 10 (roughly 10^" 
yr) starting from two initial DM profiles (see Table I and text). 
Colored curves: benchmark models; dashed lines: av = 0; 
dotted curves: av = ?, x 10"^^cm^s~^, m = 50 GeV; thick 
line: initial DM density. 

by self-annihilations but for t > lO^yr « Theat, heating 
of dark matter by stars dominates. The change in J5 
(Fig. 2b.) is dramatic, with final values in the range 
103<J<10^ 

We define the boost factor 6 as J/ Jn, with J the value 
in the evolved model and Jn the value in a p cx 
(spike-free) halo with the same density normalization at 
r — Rq. Figure 2 shows boost factors at r = 10 for each 
of the models in Table I. We found that the dependence 
of B = logj^Q b on av/m could be very well approximated 
by the function 

B{X) - B„ 



-(1/2)(B„.„ 
with X = logio(crt;/10~30 



S™„){l+tanh[C7i(X-C2)]} 
(6) 

loyuv/.^ ...xx3s-i)/(TO/100GeV). Table 
II gives values of the fitting parameters at t ~ 10 in each 
of the models with a spike. The boost factor is indepen- 
dent ofav for low av, since annihilations are unimportant 
in this limit, and also for high av, since annihilations ef- 
fectively destroy the spike. 

We now apply these results to the study of high-energy 
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FIG. 2: (a) Evolution of the dark matter density at a ra- 
dius of 10"^rfi « 2 X 10"^pc in a p ~ r~^ '^^ spike, for 
m = 200 GeV, and av = 10"^''cm^s"^ The upper(lower) 
set of curves correspond to an initial density normalization 
at Vh of 1O(1OO)M0 pc"^. In order of increasing thickness, 
the curves show the evolution of p in response to heating by 
stars; to self-annihilations; and to both processes acting to- 
gether. Time is in units of Thoat defined in the text; r = 10 
corresponds roug hly to 10^° yr. (b) Evolution of J averaged 
over an angular window of 10~^ sr. 



TABLE II: Parameters in the fitting function for the boost. 
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6.29 
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2.35 
0.38 
7.36 
3.00 
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0.41 
0.72 
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-0.15 



7-rays from dark matter annihilations at the GC. An 
early detection by the EGRET collaboration, of a 7- 
ray source coincident with the position of the SBH 41| , 
has not been confirmed by a subsequent analysis 42| . 
However, Air Cherenkov Telescopes like HESS, CANGA- 
ROO and VERITAS have all detected a source coincident 
within their angular resolution with the GC SBH. In par- 
ticular the HESS data suggest a spectrum extending up 
to 10 TeV, with no apparent cut-off It is difficult 

to interpret the observed emission as due to DM anni- 
hilation, since usual DM candidates are lighter than the 
required 10 TeV, and since the spectrum is quite fiat. 
The latter problem can be solved by considering pro- 
cesses like XX ~^ ^^7 113 1 where £ is a charged lepton, a 
channel heavily suppressed for neutralinos, but open for 
Kaluza-Klein particles. Although the contribution of the 
total fiux is small (the channel is suppressed by a factor 
a/n with respect to the annihilation to charged leptons), 
the corresponding photon spectrum is very flat, with a 
sharp cut-off at an energy corresponding to the particle 
mass. The other problem, i.e. the high dark matter par- 
ticle mass required to reproduce the HESS data, can be 
solved in the framework of some specific theoretical sce- 
narios, such as those proposed in Refs. [s^ liof l. In this 
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[4] L. Bergstrom, P. UUio and J. H. Buckley, Astropart. 
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FIG. 3: Boost factors 65 (AQ — 10^'') as a function of av/m 
at T = 10 for the DM models of Table I. Hatched region is 
the approximate boost factor required to explain the HESS 
7-ray detection if the particle mass is ~ 10 TeV and av = 
3 X 10"^'^cm^s"^ (vertical line). 



boost factor of order 10'^ ^ b < 10^ is required to 
match the observed normaUzation, for particle masses of 
order 10 TeV and cross sections of order 3 x 10~^^cm'^s~^. 
Figure 2 shows that such boost factors are achievable in 
the adiabatically compressed DM models, p ^ r~^-^, es- 
pecially if a spike is initially present, although the spike is 
not required. We note that the particle models discussed 
above could easily evade the synchrotron constraints dis- 
cussed in_0, [ij, Looking for example at Fig. 6 
of Ref. [46j, we note that the synchrotron constraint is 
weaker for heavier masses, and the annihilation rate, in 
the case of the evolved Agp profile discussed above, is 



suppressed by many orders of magnitude with respect to 
the case discussed in |46j| , corresponding to a non-evolved 
Nsp profile (Table I). A detailed analyis of indirect de- 
tection of supersymmetric and Kaluza-Klein DM in light 
of this work will be presented elsewhere. 
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